Migration and proliferation dichotomy in tumor cell invasion 



Sergei Fedotov 1 and Alexander Iomin 2 
1 Department of Mathematics, University of Manchester, Manchester M60 1QD, UK 
2 Department of Physics, Technion, Haifa, 32000, Israel 

We propose a two-component reaction-transport model for the migration-proliferation dichotomy 
in the spreading of tumor cells. By using a continuous time random walk (CTRW) we formulate 
a system of the balance equations for the cancer cells of two phenotypes with random switching 
between cell proliferation and migration. The transport process is formulated in terms of the 
CTRW with an arbitrary waiting time distribution law. Proliferation is modeled by a standard 

\^ • logistic growth. We apply hyperbolic scaling and Hamilton- Jacobi formalism to determine the 

overall rate of tumor cell invasion. In particular, we take into account both normal diffusion and 
anomalous transport (subdiffusion) in order to show that the standard diffusion approximation for 

O^l ■ migration leads to overestimation of the overall cancer spreading rate. 
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PACS numbers: 05.40.-a, 05.40.Fb, 87.15. Vv, 87.17.Ee, 82.39.Rt 
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C ' Extensive investigations have been devoted to the modeling of cancerous growth (see, for example, reviews [USEl- 
and references therein). Although a great deal of progress has been made in this theory, especially for solid tumors, 
for which growth is basically due to cell proliferation, our understanding of malignant gliomas, the diffusive and highly 
invasive brain tumors, is much less complete (see, a review 0). The main reason for this is that unlike solid tumors, 
gliomas not only are able to proliferate but also to invade the surrounding brain parenchyma actively. The surgical 
resection of diffusive tumors is ineffective since the cancer cells have already invaded the surrounding brain tissue. 
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This leads to recurrence of tumor, and the prognosis for patients suffering from malignant gliomas is very poor. Thus 
proliferation and especially migration of gliomas provide a significant challenge for modelling, and this is why the 
invasiveness of tumors has been studied extensively in recent years (see, for example, 0,0,0)- 
i ^ , '■ Invasion, itself, is a very complex process of receptor- mediated transport 0, which involves several steps of cell 
migration and proliferation (see a review 0). Experimental evidence indicates the lower proliferation rate of migratory 
OnJ i cells in comparison with the tumor core, which indicates an inverse correlation between mobility and proliferation of 
cell population. The existence of this important phenomenon was supported by numerous experimental data obtained 
in vitro and clinical data obtained in vivo A\. It was formulated by Giese et al. as a migration-proliferation 
dichotomy. It turns out that proliferation and migration of tumor cells are mutually exclusive phenotypes: the 
spreading suppresses cell proliferation and visa versa. The molecular mechanism for this dichotomy has been suggested 
' m ' an d then an active implementation for the numerical modelling of the brain tumor and its fractional topology 
\Q \ has been established Q. It turns out that this behavior of cells is an inherent process of a so-called continuous time 
• random walk (CTRW). This transport concept, based on j ump and waiting time distributions, has been extensively 
and successfully employed for numerous applications [l(J,[ll|]. Migration-proliferation dichotomy was formulated in the 
framework of the CTRW in 12] . The primary focus was on the influence of cell fission on transport properties of cells. 
An essential decrease in cell motility during fission time, or their self-entrapping, is determined by the interaction of 
cells with their environment. In vitro experimental observations of cell transport confirm the essential decrease in cell 
motility during cell proliferation [l3j . 
. £h ! Usually the random mobility of tumor cells is described by Fick's law. However, it has been shown that the 
diffusion approximation for the transport process together with a logistic growth yields an overestimation of the 
overall propagation rate |l4L ll5| . Since the tumor cells' migration is the most critical feature of brain cancer, causing 
treatment failure, the transport has to be properly understood. Therefore we need to extend the diffusion analysis by 
introducing a more realistic description of the transport of mobile tumor cells. It is one of the main purposes of this 
Letter to take into account anomalous transport (subdiffusion) leading to slow mobility of cancer cells in the invasive 
zone. 

In this Letter we propose an alternative approach for the migration-proliferation dichotomy. We employ a two- 
component CTRW, assuming that the glioma cells are of two phenotypes. In state 1 (migratory phenotype) the cells 
randomly move but there is no cell fission. In state 2 (proliferating phenotype) the cancer cells do not migrate and only 
proliferation takes place. The exact mechanism of switching between the two phenotypes is not known. An interesting 
deterministic mechanism for this phenotype switch has been suggested recently in |l6|. However, the mathematical 
modeling involves many parameters, some of which are difficult to estimate. Here we propose the stochastic approach 
for the proliferation-migration switching that involves only two parameters. We assume that the cell of type I remains 
in a state I during a waiting time t% and then switches to a cell of type 2. After a waiting time T2, spent in a state 2, 
it switches back to a cell of type f . Both waiting times n and Ti are mutually independent random variables. In this 
paper we consider them exponentially distributed with parameters fa and fa: namely, /(r^) = fa exp (— fa^i). Here 
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the parameters /3, are the switching rates, namely, 0\ is the switching rate from state 1 to 2, while 2 determines the 
transition rate 2 — > 1. The last parameter can control the lower migratory cell proliferation. Indeed, if the average 
time (t 2 ) — 1//32 is much less than time of cell doubling then the proliferating cells in state 2 do not have enough 
time to proliferate even with the same high rate as the core of cells. 

Let us introduce the density for the cells of migratory phenotype (cells of type 1), Ui(t, x), and for the cells of 
proliferating phenotype (cells of type 2), n 2 (t,x). The balance equations can be written as follows 



ni(t,x) = ni(0,x)^(t)e 0lt + J J m(t - s,x - z)p(z)tp(s)e 0lS dzds 

+02 / n 2 (t- s,x)*(s)e-' 9lS ds, (1) 
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(t,x) = n 2 {0,x)e' l32t + U n 2 (t - s, x)(l - n 2 (t - s, x)/K)e- 02S ds 



+0i I ni{t - s,x)e~ 02S ds, (2) 
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where p(z) is the probability density for migration jump length, while ip(s) is the probability density of waiting times 
between jumps, and ^(t) = 1 — f ip(s)ds is the probability that a cell of type 1 makes no jump until time t. The 
exponential factor e~ 0it = J°° /(r,)dr, is the probability that cells of phenotype i do not switch until time t. Eq. (JTJ 
describes the balance of cells of type 1 at time t at position x. The first term on the right hand side of the equation 
represents those cells of type 1 that stay up to time t at position x such that no jump occurs, and no switch 1 — ► 2 
takes place. The independence of the random jumps and switching gives us the probability ^(^e - ^ 1 * while the first 
factor ni(0, x) is the initial density of cells of type 1. The second term describes the number of cells of type 1 arriving 
at x up to time t due to the following random mechanism of migration: the cell of type 1 at time t — s at position 
x — z waits a random time s before jumping at a distance z and remains a cell of type 1. This process is determined 
by the transition probability ip(s)p(z). The limits of the space integral are determined by the boundaries. The last 
term in Eq. Q represents the number of cells of type 2 that switches to the cell of type 1 up to time t and remaining 
the cells of type 1 (due to the factor e~ 0lS ). It also takes into account the fact that if transition 2 — > 1 happens at 
time t — s, then no jump takes place during the remaining time s (due to the factor 

Regarding Eq. J5J), the first term on the right hand side has the same physical meaning as one in Eq. QJ. The 
second term is the logistic growth 01 f° r cens 01 type 2, which occurs providing that no switch takes place up to time 
t. Here U is the cell proliferation rate and K is the carrying capacity of the environment. The last term of Eq. 
represents the number of cells of type 1 switching to the state 2 over the time interval (0, t). Note that one-component 
balance equation involving transport and production term has been analyzed in |15| . 

The balance equations JIJ and J5J can be written as the system of integro-differential equations. By using the Laplace 
transform hi(H) — J °° e~ Ht rn(t)dt and presenting the left hand side of the equations in the form Hhi{H) — hi(0) 
which is the Laplace transform of the time derivative, one obtains 

—q£ = J oc(t — s) J [rii(s,x — z) — rii(s,x)] p(z)dzds — 0ini + 2 n 2 , (3) 

^ =Un 2 (l-n 2 /K) + ir ii- 2 n 2 , (4) 
at 

where the 'memory' kernel a(t) is defined in terms of its Laplace transform 

a(H) = (JL±MM±M , (5) 



(l-^iT + A)) 



with tf){H) = f °° ip(t)e~ Ht dt. Note that the eq uivalence of one-component balance equation to a master equation 
involving memory kernel has been shown in [lOj. It should be emphasized that it is impossible to find an explicit 
expression for memory kernel a(t) for arbitrary choices of waiting-time pdf tp(t). In what follows we will be concerned 
with the overall rate of the spreading of gliomas. It turns out that this rate depends on the Laplace transform a(H) 
rather than a(t). That is why the formula iJSJ plays a crucial role in this Letter. 
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It is natural to assume that jumps of migrating cells are small and there is no convection (J zp{z)dz = 0). Expanding 
rii(s, x — z) in the Taylor series in Eq. Q, we obtain 



dni 



d 2 ni 



(6) 



where a 2 — J z 2 p(z)dz. Generalization on 3D is straightforward, namely, the second derivative is replaced by the 
Laplace operator A. 

Now we are in a position to find the overall rate at which the cancer cells spread. The main purpose here is to find 
the dependence of the rate of invasion on the statistical characteristics of the random switching process, fi\ and /3 2 , 
and random walk in space, a 2 and ip(t). We expect that the system of equations together with appropriate initial 
conditions has a traveling wave solution (planar front) with some velocity u common to both densities n\ and n 2 . The 
objective here is to find the rate u without resolving the shape of the traveling waves |l4lll8j |. For this purpose we use 
a hyperbolic scaling x — > x/e, t — > t/e and the rescaled densities nf (t, x) — rii (t/e,x/e). We apply the exponential 
transformation 



n\ (t, x) = Ai exp 



G (t, x) 



1,2. 



(7) 



where positive constant A\ and A 2 represent the stable equilibrium points of the densities n\ and n|. Our purpose 
is to find an equation for G (t, x) which gives us the spreading front position x (t) in the limit of the long-time and 
large-distance, from the equation G(t,x(t)) = 14]. To ensure the minimal spreading rate we use the front-like 
initial conditions: rij(0, x) = Ai for x < 0, and rij(0, x) — for x > 18]. Substituting Q into the equations for 
the densities n\ and n§, one obtains two equations for A\ and Ai in the limit e — > 0. This system has a non-trivial 
solution when the corresponding determinant is equal to zero. This yields a generalized Hamilton- Jacobi equation, 
involving two first derivatives dG/dt and dG/dx : 



a 2 fdG 



l ~ 1+ yU 




%t s i)(s)e-P ls ds 



1-U e^ s e^ 8 ds 



Pilh / e^ s 4'(s)e-' 3lS ds x / e^ s e-^ s ds = 



(8) 

Note that inferring Eq. |8|. we do not make any assumptions regarding waiting time pdf ip(t). If we introduce the 
Hamiltonian function H = —dG/dt, the generalized momentum p = dG/dx, and the Laplace transform ip(H) = 
J °° ip(t)e~ Ht dt, then the Hamilton- Jacobi equation JHJ takes the form 
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/?i/3 2 (i-y>(g + /3i)) 



The latter equation is important, since it allows us to find the overall spreading rate u 



H 

pJh)- 



dp 
~d~H 



P(H) 
H 



(9) 



mintf {H/p(H)} by using 



(10) 



In the symmetrical 31? case, Eq. (J^J corresponds to the Hamiltonian motion in the radial direction. Let us illustrate 
the use of the above theory through two typical distributions for the waiting-time pdf ip(t). 

First, we consider a probability distribution function for the exponentially distributed waiting times: = 
T _1 e _t / T . We find 4>{H) = (1 + Ht)^ 1 and a{H) = r _1 , and therefore a(t) = r~ 1 5(t). This corresponds to 
the classical Fick's law for transport with the diffusion coefficient D = a 2 /2r. Thus we have a classical system of 
reaction-diffusion equations such that the equation for the migratory cells is 



dn\ 
~~dt 

The momentum p(H) can be found from 



D 



d 2 ni 
dx 2 



D 



Pini + /3 2 n 2 . 



Pith 
D(H + 02-Uy 



(11) 



(12) 
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If we assume that Pi — P2, we can find from (|10|) and 1|12[1 p = (U/D) 1 ^ 2 , and H = U. Therefore, the spreading rate 

1/2 

is Mo = (UD) which is half of the classical Fisher-KPP propagation speed. This is a very interesting result showing 
that the propagation rate is independent of the random migration-proliferation switching when cell transport is the 

1 /2 

Brownian motion and Pi = (32- When Pi ^ 02 one can find the ratio of the propagation rate u and uq = (UD) as 

H 2 (H + (3 2 -U) 



u J U ((H + 02 — U)(H + Pi) - /3iP 2 ) ' 

The situation changes for the power law distribution (anomalous transport): ip(t) ~ (r/i) 1+7 with < 7 < 1. 
This is the case when the mean waiting time is divergent: < t >= 00. This assumption alone leads to the temporal 
fractional differential operator and corresponding anomalous diffusion equation . The mean squared displacement 
for mobile cells is 

where Z? 7 = ct 2 /2t 7 is the generalized diffusion coefficient with the dimension cm 2 s" 7 . One of the main aims of this 
Letter is to find the overall propagation of cancer cells as a result of interaction of the anomalous migration 1|14[) . 
logistic proliferation and random migration-proliferation switching. For this purpose it is more convenient to define 
if>(t) by its Laplace transform 4>{H) = (1 + (Ht) 1 ) 1 [nj, such that the momentum p(H) can be found from @ 

2 = (g + Ap _ PifajH + ^y- 1 

P L» 7 D J (H + p 2 -U)' 1 ' 

This formula together with I|1U|) allows us to find the overall propagation rate of tumor cells u 7 in the fractional 
diffusion case. It is clear that the case 7 = 1 corresponds to the normal diffusion approximation for cell migration 
(see H12|)l. One can find from IjlUI) . (|12[l and (|I5|) the ratio of the anomalous propagation rate u 7 and the normal rate 
u determined by (|13|) : 



u 



(ff 7 r + Ar)^, (16) 



where H 1 is the solution of dp/dH = p(H)/H. Since the "microscopic" time r is much smaller than the "mesoscopic" 
reaction time U^ 1 and switching time /Jf 1 and i? 7 ~ U, we conclude that Ht + Pit < 1. It follows from (|16|l that 
the ratio w 7 /w increases with 7 in the interval < 7 < 1. This means that the standard diffusion approximation 
leads to overestimation of the overall cancer spreading. It is clear from these two examples of normal and anomalous 
diffusions that the advantage of balance Eqs. Q and J5J is that they are related to "mesoscopic" description of 
migratory cancer cells, and give us the statistical meaning of the reaction-diffusion equations or fractional equations 
that are introduced usually phenomenologically 

In summary, we present a two-component model for a migration-proliferation dichotomy in the spreading of tumor 
cells in the invasive zone. We use a probabilistic approach based on the CTRW theory for migration, logistic growth 
and random proliferation-migration switching with exponentially distributed waiting times. Our approach is not 
restricted to the specific mechanism of proliferation describe d by a logistic growth. Moreover, Eq. © for proliferation 
can be accompanied by a nutrient control or chemotaxis |19j . The main point of the paper is that cancer cell 
transport is subdiffusive rather than diffusive described by Fick's law (the cancer cells are not Brownian particles!). 
The advantage of our approach is that it allows us to take into account anomalous (subdiffusive) transport within the 
general scheme of migration, proliferation and phenotype switching. We show the equivalence of balance equations to 
a system of master equations involving memory kernels for the transport of mobile cells. By using a hyperbolic scaling 
and Hamilton- Jacobi formalism we derive formulae for the overall spreading rate of cancer cells. We show that the 
memory effects (subdiffusion) leads to a decrease in propagation rate compared to a standard diffusion approximation 
for transport. An analytical expression for the memory kernel can be obtained for more complicated processes. For 
example, for the family of gamma distributions with parameters m and r, ip(t) = T -' m t m - 1 e - t / T /T(m) [2(J]. We have 
$[H) = (1 + H T y m and 

~ (m (H + Pi) 



(1+Ht + p x r) m - 1 



For example, if m = 2, then a(H) = r 1 (2 + Ht + Pit) 1 and the memory kernel is ait) = r 2 e ( 2 +i 3 i T ) t / T . The 
integro-differential Eq. © can be written as the hyperbolic reaction-transport equation, and corresponding traveling 
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wave solutions can be found in 2l| . Renovation processes with arbitrary probability densities for switching waiting 
times will be considered in the future publications. 



[1] N. Bellomo and L. Preziosi, Math. Comp. Model. 32, 413 (2000). 
[2] K. Swanson et al. J. Neur. Sci. 216 1 (2003). 

[3] Biomedical Applications ed. by H. Byrne, in Modeling and Simulation in Science, Engineering and Technology (Birkhauser 

Boston and Basel, 2006), in press. 
[4] A. Giese et al. J. Clin. Oncology 21 1624 (2003). 

[5] L. M. Sander, Th. S. Deisboeck, Phys. Rev. E 66, 051901 (2002); S. Habib, C. Molina-Paris, Th. S. Deisboeck, Physica A 
327 501 (2003). 

[6] T. Bollenbach, K. Kruse, P. Pantazis, M. Gonzales-Gaitan, and F. Jiilicher, Phys. Rev. Lett. 94, 018103 (2005). 

[7] A. Giese et al, Int. J. Cancer 67, 275 (1996). 

[8] A. Wells, Int. J. Biochem. Cell Biol. 31, 637 (1999). 

[9] Y. Mansury and T. Deisboeck, Physica A 331, 219 (2004); Physica D 196, 193 (2004). 
[10] E. W. Montroll and B. J. West, in Fluctuation Phenomena, edited by E. W. Montroll and J. L. Lebowitz (North-Holland, 
Amsterdam, 1979); W. Montroll and M. Shlesinger, On the Wonderful World of Random Walks (Elsevier Science Publishers 
BV, 1984). 

[11] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000). 

[12] A. Iomin, J. Phys.: Conference Series 7, 57 (2005); A. Iomin, WSEAS Trans. Biol. Biomed. 2, 82 (2005); A. Iomin, Phys. 

Rev. E 73, 061918 (2006). 
[13] P. Dieterich, private communication. 
[14] S. Fedotov, Phys. Rev. Lett., 86, 926 (2001); 
[15] S. Fedotov, and V. Mendez, 66, 030102 (2002). 

[16] C. A. Athale, Y. Mansury, and T. S. Deisboeck, J. Theor. Biol. 233, 469 (2005). 
[17] J. D. Murray, Mathematical Biology (Springer- Verlag, Berlin, 1989). 

[18] M. Freidlin, Markov Processes and Differential Equations: Asymptotic Problems (Birkhauser, Basel, 1996). 
[19] E. Khain, L. D. Sander, and A.M. Stein, Complexity 11, 53 (2005); E. Khain, L. D. Sander, Phys. Rev. Lett. 96, 188103 
(2006). 

[20] S. Fedotov, and Y. Okuda, Phys. Rev. E. 66, 021113 (2002). 
[21] J. Fort, and V. Mendez, Rep. Progr. Phys. 65, 895 (2002) 



